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ABSTRACT 

We present a power spectrum analysis of the Herschel-SPWJL observations of the Polaris flare, a high Galactic latitude cirrus cloud midway 
between the diffuse and molecular phases. The SPIRE images of the Polaris flare reveal for the first time the structure of the diffuse interstellar 
medium down to 0.01 parsec over a 10 square degrees region. These exceptional observations highlight the highly filamentary and clumpy structure 
of the interstellar medium even in diffuse regions of the map. The power spectrum analysis shows that the structure of the interstellar medium 
is well described by a single power law with an exponent of -2.7 ±0.1 at all scales from 30" to 8°. That the power spectrum slope of the dust 
emission is constant down to the SPIRE angular resolution is an indication that the inertial range of turbulence extends down to the 0.01 pc scale. 
The power spectrum analysis also allows the identification of a Poissonian component at sub-arcminute scales in agreement with predictions of the 
cosmic infrared background level at SPIRE wavelengths. Finally, the comparison of the SPIRE and IRAS 100/im data of the Polaris flare clearly 
assesses the capability of SPIRE in maping diffuse emission over large areas. 
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1. Introduction 

Understanding the physical mechanisms involved in the forma- 
tion of stars is still one of the major aspirations of modern as- 
trophysics. One of the many aspects of this problem is related 
to the organisation of matter in the diffuse interstellar medium 
(ISM) as its structure provides the initial conditions for the for- 
mation of molecular clouds. The physics that impacts on the 
structure of the ISM is complex : it involves compressible and 
magnetic turbulent motions in a thermally unstable flow. Even 
though turbulence organises the flow up to scales of tens of par- 
sec ( Elmegreen & Scalo 2004 1 , important non-linear physical 
processes take place at the sub-parsec scales. Being able to study 
the kinematics of the gas at such scale is essential but the density 
structure is also a key tracer of the physics involved. 

Dust emission from uniformely heated regions of the local 
interstellar medium (i.e. exempt of local production of UV pho- 
tons by stars) allows us to directly trace the density structure 
of the ISM. With its unprecedent angular resolution, the ESA 
Herschel Space Observatory (Pilbratt et al. |2010 l allows us for 
the first time to study the structure of the ISM at the 0.01 pc scale 
over several square degrees. We present a first analysis of the 
250, 350 and 500 //m emission maps of the Polaris flare obtained 



with the SPIRE instrument (| Griflin et al.||2010l l. These observa- 
tions are part of the "Evolution of Interstellar dust" ( Abergel et 
aTlOTO) and "Gould Belt survey" ( [Andre et al.||2010[ ) SPIRE 
key projects. 

The Polaris flare is a high Galactic latitude cirrus cloud, 
located at a distance < 150 pc (|FaIgarone et al.| I998|l. It 



has significant CO emission ([Heithausen & Thaddeus|" 1990 



Meyerdierks & Heithausen 1996 Falgarone et al. 1998 1 and 
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colder dust grains than typical diffuse clouds ( Lagache et al. 
1998t [Bernard et al.|ri999| l. As the Polaris flare does not show 
any sign of star-formation activity (only pre-stellar cores are de- 
tected with SPIRE - [Andre et ah] ( |20I0[ l, it is the archetype of 
the initial phases of molecular cloud formation. We present a 
power spectrum analysis of the dust continuum emission from 
big grains using SPIRE observations. These data allow us to 
study homogeneously the structure of the diffuse ISM from 0.01 
to 8 pc. 



2. SPIRE observations of the Polaris flare 

The 250 fim SPIRE map, covering 10 square degrees of the 
brightest part of the Polaris flare, is shown in Fig. [T] The field 
was observed twice in nearly peipendical directions to optimize 
the restoration of diffuse emission. The SPIRE data were pro- 
cessed with HIPE (version 2.0) applying standard corrections 
for instrumental effects and glitches. The 1/f noise component 
was removed using the "temperature drift correction" module, 
and naive maps were computed. 
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Fig. 1. SPIRE 250 fim map of the Polaris flare. Units are in MJy sr ' . The zero level was set by correlation with the IRAS/IRIS 
100 fim data. 



This map reveals for the first time the structure of the dif- 
fuse interstellar medium on scales ranging from 0.01 to 8 pc. 
Compared to the previous vision of the diffuse interstellar 
medium given by lower resolution observations (e.g. IRAS) 
these observations reveal a structure with strong contrast at small 
scales. Numerous small scale clumps are seen in the map even 
in the most diffuse regions (see examples in the Appendix). This 
high-resolution view of the diffuse ISM also reveals its highly fil- 
amentary structure with narrow threads of matter following the 
larger scale organisation. These observations bring new insights 
into the small scale structure of the ISM, and they will certainly 
help understand the physical processes dominating the dynami- 
cal evolution of matter towards the formation of stars. This task 
is obviously beyond the scope of the present paper. 



3. Comparison with IRAS 100 jum: cliecl< of diffuse 
emission restoration and dust spectrum 

The estimate of the power spectrum of the interstellar medium 
emission can only be done with observations that restore the 



power observed at all scales. This implies a great control of in- 
strumental effects that could affect the baseline (additive effect) 
or the gain (multiplicative effect) of the detectors over the whole 
period of the observations. In order to assess the quality of the 
diffuse emission restoration by SPIRE over such a large field 
we made a comparison over the whole field with the 100 /urn 
IRAS (IRIS) data ( [Miville-Deschenes & Lagachej [2005] l which 
are known to have a good description of the interstellar emission 
at all scaleqj The main limitation of this exercise is the differ- 
ence in wavelength between SPIRE and IRAS, but even though 
local variations of the dust emission spectrum are expected, the 
fact that both datasets are dominated by the emission from the 
big grain population is instructive. 

We performed the following linear regression fit : S(A) = 
G X 5(100) + So, where 5(100) is the 100 jum IRAS/IRIS map 
from which the average value of the cosmic infrared background 
at 100 yum (0.7 MJy/sr -[Miville-Deschenes et aL](|2007[)) was re- 



' in IRIS the variation with scale of the IRAS detectors gain was 
corrected and the emission at scales larger than 30' was made consistent 
with DIRBE, which was designed for full control on systematic efi'ects. 
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moved, and S{A) is the SPIRE map at wavelength A, convolved 
to the IRAS resolution (4.3 arcmin) and projected onto the na- 
tive 1.5' grid of IRAS. The regression coefficients G and offsets 
^o found at each SPIRE wavelength are given in Table [T] Even 
though the correlations are good (correlation of 0.85), there is 
significant variation around the linear fit. Looking at the dif- 
ference map (i.e. S{A) - G x 5(100)) localized variations are 
seen, which reflect expected modifications of the dust emission 
spectrum. These small and intermediate scale variations sit on 
a fainter large scale structure uncorrected with the emission 
and therefore probably unrelated to variations of dust properties. 
This residual, which is less than 10% of the large scale emission 
fluctuations, could be attributed to residual imperfections in the 
data processing. 

The factors G obtained for the SPIRE-IRAS correlation can 
be used to estimate the average dust emission spectrum in the 
field. The fit of a grey bodjj^to the IRAS-SPIRE correlation co- 
efficient shown in Fig.[2]gives Tj = 14.5 ±1.6 and/? - 2.3 +0.6, 
in agreement with what was measured by Bernard et al. ( I999| l 
using PRONAOS and ISOPHOT data on a 30' x 6' region in 
the brightest part of the Polaris flare {Tj = 13.0 + 0.8 and 
j6 = 2.2 + 0.3). This provides a first order sanity check of the 
quality of the SPIRE gain calibration. In addition, looking at 
(S(A) - So)/S(lQQ), we find no systematic correlation of G with 
intensity which agrees with the fact that the SPIRE diffuse emis- 
sion calibration is not scale-dependent at scales larger than the 
IRAS beam. 



4. Power spectrum analysis 

The power spectrum of the 250 yum map of the Polaris flare, con- 
verted to Jy^/sr is shown in Fig. [3] The black dots in the bottom 
plot show the power spectrum computed on a 2.8° x 3.1° area of 
the map where all data points are defined and from which bright 
point sources were removed. An apodization factor of 0.97 was 



applied prior to the Fourier Transform ( Miville-Deschenes et al. 
[2002) 1. 

The power spectrum is typical of infrared emission of high 
Galactic latitude fields with a power-law type spectrum con- 
volved by the instrument transfer function ((p), and a flat noise 
part (A^) at high k. The power spectrum is modeled accordingly : 
P(k) = (p{k)Psky(k) + N(k). The white noise term stands out very 
clearly in all power spectra. Its level is estimated as the average 
of P{k) for 0.15kmax < k < k^ax, where k^ax is the maximum k 
available (i.e. twice the pixel size - see TablefTb. At each SPIRE 
wavelength this corresponds to scales smaller than the beam size 
where the noise dominates. The recovered noise levels are given 
in Tablefl] They are comparable with the expected sensitivity for 
two repeats. 

The dark blue dots in Fig. l3]show the power spectrum of the 
250 yum map, noise removed and divided by the transfer function 
if>{k) estimated using the official SPIRE beam profiles obtained 
on observations of Neptune. We emphasize here that the SPIRE 
beam shapes cannot be approximated by a Gaussian for the level 
of precision needed in this power spectrum analysis. Not taking 
into account the secondary lobes of the transfer function with a 
Gaussian would produce an artificial break in the power spec- 
trum at scales of 0.1-0.2 arcmin ' with a steepening of the slope 
at small scales. 

Once corrected for noise and (p{k) the power spectrum shows 
a rather straight power law with a sUght flattening at wavenum- 
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Fig. 2. Average big dust emission spectrum in the Polaris flare 
as determined by correlation of the SPIRE and 100 fxm data over 
the whole field. In this plot the correlation coefficients (G in 
Table [Til were scaled to the average 100 /urn brightness in the 
SPIRE field. The fit of the big grain emission spectrum gives 
Td = 14.5 ± 1.6 andjS = 2.3 ± 0.6. 



bers larger than 1 arcmin ' , typical of a white component due to 
point sources and the unresolved cosmic infrared background 
(CIB) Poissonian fluctuations. To extract the interstellar con- 
tribution to the power spectrum we fitted ffie dark blue curve 
with Psky(k) = Ajsujk^ + Pq, where Pq is the white level due to 
point sources and the CIB. The fit was done on scales between 
k=0.025 arcmin"', to exclude the largest scales where the IRIS- 
SPIRE comparison showed significant differences, and about 
twice the FWHM (2.0, 1.2 and 0.86 arcmin"' at 250, 350 and 
500 fim respectively) to exclude scales contaminated by residual 
noise. The recovered y and Pq are given in Table [T] The uncer- 
tainties quoted in Tablefl] include both the statistical error of the 
fit and a bias of a ±10% uncertainty on A^. The increase of the 
uncertainty for y from 250 to 500 yum is due to the range of scales 
on which the fit is done, which get shorter with increasing beam 
size. Even though the uncertainties for Pq are significant, it is 
remarkable that the levels obtained agree very well with the pre- 
dictions of |Fernandez-Conde et al.| ( |2008| ) for the CIB fluctuation 
level at the SPIRE wavelengths. 

At 250 fim the slope of the power spectrum is y = -2.65 ± 
0.10 at wavenumbers from A:=0.025 to 2.0 arcmin"'. In the top 
panel of Fig.[3]the green dots represent the power spectrum (mul- 
tiplied by k^ ) of the SPIRE 250 jim map corrected for A^, </> and 
Pq, leaving only the interstellar contribution. Within the uncer- 
tainties, the slopes are similar at all wavelengths (see Table [TJ. 
They are also compatible with the IRAS/IRIS 100 fim power 
spectrum measured on 12° x 12° region centered on the SPIRE 
fielqj(red dots). However we note that the SPIRE power spec- 
trum is systematically lower than IRAS for the first few large 
scale modes, in agreement with what was shown in Sect.l3] This 
effect, which could be partly due to suppresion of large scale 
modes in the SPIRE map-making, certainly deserves further in- 
vestigation especially by comparing data in the common bands 
of SPIRE and Planck. 



5. Discussion and conclusions 

The power spectrum analysis of the 250, 350 and 500 jum SPIRE 
maps of the Polaris flare described here shows ffiat ffie -2.7 slope 



measured at large scales using IRAS (Miville-Deschenes et al. 



^ ly on By(Tj)v ^, where T,/ and fj are the big grain temperature and 
emissivity index, and S,, is the Planck function. 



^ The IRAS/IRIS power spectrum shown here has been corrected for 
noise, the transfer function and the CIB fluctuation level. It was scaled 
to match the SPIRE power spectrum at 10-15 arcmin scales. 
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1.7±0.3 (0.008) 
0.7±0.1 (0.003) 


20.1±3.1 (0.1) 
10.1±1.6(0.05) 
4.3±0.7 (0.02) 


1.26 
0.55 
0.34 


-2.65 ±0.10 
-2.69 ±0.13 
-2.62 ±0.17 


5±2 X 10' 
4±2 X 10^ 
1±1 X 10^ 



Table 1. SPIRE observations of the Polaris flare. Column 4 and 5: gain and offset coefficients of the SPIRE-IRIS 100 jum 

correlations. The uncertainty on G represents the rms of the ratio S (A) / S (lOOfim) once the offset So is removed from the SPIRE 
data. Similarly the uncertainty on So is the rms of S(A) - G x S(lOOfim). These two uncertainties are correlated, but they give 
more realistic estimates compared to the statistical ones obtained with the linear regression fit (given in brackets). Column 6: noise 
level estimated directly on the power spectrum for k > Q.15k,„ax (3.75 < A: < 5 arcmin"' at 250 yum). Column 7: power spectrum 
exponent of the interstellar component. Column 8: level of the Poissonian component associated with point sources and the CIB 
fluctuations. The fit of the power spectrum (y and Pq) was done between k=0.025 arcmin"' and about twice the FWHM (2.0, 1.2 
and 0.86 arcmin' at 250, 350 and 500 /im respectively). 
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Fig. 3. Power spectrum of the SPIRE 250 jjm map of the Polaris 
flare. Bottom: The black dots are the raw power spectrum (com- 
puted with an apodization of 0.97). The horizontal black line is 
the white noise level estimate. Blue dots are the power spectrum 
noise removed and divided by the psf estimated on Neptune ob- 
servations. The green horizontal line and the dashed line are the 
source and interstellar components estimated from the blue dots 
power spectrum respectively, on scales 0.025 < k <2 arcmin-'. 
Top: The green dots are the power spectrum of the interstellar 
component (i.e. corrected for noise, psf and sources) multiplied 
by k^-^^. The red dots are the power specti'um of the IRAS/IRIS 
100 fim emission in a 12° x 12° region centered on the SPIRE 
field. The dashed line is the same as in the bottom figure. 



2007| l or CO data (J Stutzkiet aL]|I998[ ) extends by more than one 
order of magnitude in scales, down to 30". The high-resolution 
data obtained here on a 10 square degree field are in accordance 
with the results of [Heithausen et al.| ( |I998| l, who also found 
a slope of -2.7 on a similar range of scales combining high- 
resolution IRAM data of two 6' by 8' fields with lower reso- 



lution CfA and KOSMA data of a larger field. The analysis also 
reveals a Poissonian component at a level that agrees with CIB 
estimates. 

As Polaris is uniformely heated by the interstellar radiation 
field, the dust grain temperature is likely to be rather uniform 
over the field. In these conditions brightness flutctuations in the 
SPIRE bands are dominated by dust column density fluctuations 
and the power spectrum of dust emission probes directly the 
power spectrum of the dust volume density in three dimensions. 
Assuming gas and dust are well mixed in this diffuse cloud, the 
power spectrum of the gas volume density would then have a 
slope of -2.7 ±0.1 down to 0.01 pc. The constant slope of the 
power spectrum suggests that the turbulent cascade is still the 
main agent organizing the structure of the ISM at the 0.01 pc 
scale. 

Finally it is important to point out that at this early phase of 
the Herschel mission it is impossible to rule out completely that 
the analysis is affected at some level by instrumental or process- 
ing effects. Nevertheless the comparison with IRAS/IRIS 100//m 
and the use of the psf measured in-flight provide a coherent pic- 
ture of the power spectrum analysis and show the capabilities of 
SPIRE at performing mapping of diffuse emission. 
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Fig. 4. Zoom on specific regions of the Polaris flare. Spire 250 fim data in MJy sr 
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Fig. 5. Zoom on specific regions of the Polaris flare. Spire 250 fim data in MJy sr 



Miville-Deschenes, M.-A. et al.: Herschel-SPIRE observations of the Polaris flare 



27.0 



T3 



26.5 



26.0 




65 



30 



125.8 125.6 125.4 125.2 125.0 124.8 124.6 

I (degree) 



Fig. 6. Zoom on specific regions of the Polaris flare. Spire 250 fim data in MJy sr 



